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The Evolution of Isotope Ratios in the Milky Way Galaxy 
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ABSTRACT 

Isotope ratios have opened a new window into the study of the details of stellar 
evolution, supernovae, and galactic chemical evolution. We present the evolution of 
the isotope ratios of elemental abundances (from C to Zn) in the solar neighbourhood, 
bulge, halo, and thick disk, using chemical evolution models with updated yields of 
Asymptotic Giant Branch (AGB) stars and core-collapse supernovae. The evolutionary 
history of each element is different owing to the effects of the initial progenitor mass 
and metallicity on element production. In the bulge and thick disk the star formation 
timescale is shorter than in the solar neighbourhood, leading to higher [a/Pe] ratios. 
Likewise, the smaller contribution from Type la supernovae in these regions leads to 
lower [Mn/Fe] ratios. Also in the bulge, the abundances of [(Na, Al, P, CI, K, Sc, Cu, 
Zn)/Fe] are higher because of the effect of metallicity on element production from 
core-collapse supernovae. According to our predictions, it is possible to find metal- 
rich stars ([Fe/H] > — 1) that formed in the early Universe as a result of rapid star 
formation. The chemical enrichment timescale of the halo is longer than in the solar 
neighbourhood, and consequently the ratios of [(C, F)/Fe] and 12 C/ 13 C are higher 
owing to a significant contribution from low- mass AGB stars. While the [ev/Fe] and 
[Mn/Fe] ratios are the same as in the solar neighbourhood, the [(Na, Al, P, CI, K, 
Sc, Cu, Zn)/Fe] ratios are predicted to be lower. Furthermore, we predict that isotope 
ratios such as Mg/ 25,26 Mg are larger because of the contribution from low-metallicity 
supernovae. Using isotopic ratios it is possible to select stars that formed in a system 
with a low chemical enrichment efficiency such as the satellite galaxies that were 
accreted onto our own Milky Way Galaxy. 

Key words: Galaxy: abundances — Galaxy: evolution — stars: abundances — stars: 
AGB and post-AGB — stars: supernovae 



1 INTRODUCTION 

Elemental and isotopic abundances are the fossils of galactic 
archaeology. Different elements are produced from stars on 
different timescales, therefore elemental and isotopic abun- 
dance ratios provide independent information on the "age" 
of a system and can be used as a form of "cosmic clock". 
The formation and evolutionary history of galaxies can be 
constrained in theoretical models by using the information 
cont ained in the elemental abundances observ ed in stars 
(e.g.. lTinslevlll980l : iPagell 19971 ; lMatteuccill200ll b The space 
astrometry missions (e.g., GAIA) and large-scale surveys 
(e.g., the high-resolution multi-object spectrograph HER- 
MES on the Anglo-Australian Telescope) will produce un- 
precedented information on the chemody namical structure 
of the Milky Way Galaxy. Theoretically iKobavashi et all 
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(2006, hereafter K06) succeeded in reproducing the aver- 
age evolution of major elements (except for Ti) in the solar 
neigh bourh ood (see also the recent study bv lRomano et ail 
l20ld b and IKobavashi fc Nakasatd j|2010l hereafter KN11) 
predicted the frequency distribution of elements in the Milky 
Way Galaxy depending on the location. 

Isotope ratios of elemental abundances can provide 
more constraints not only on galaxy formation and evolu- 
tion, but also on the detailed physics of Asymptotic Gi- 
ant Branch (AGB) stars and supernovae. The isotopic ra- 
tios of C, N, O, and Si have been measured in ancient 
meteorites and provide information on condit ions in the 
proto-solar nebula (e.g., lAnders fc Zinnei] [l993) . as well as 
the composit ion of stellar winds from pre-solar grains (e.g., 
Zinncr 1998). Attempts have been made to use meteoritic 
data to trace the chemical enrichment history of the Milky 
Way, although this has been do ne only for a limited num- 
ber o f elements including O, Mg dNittler et al.1120081 ; iNittlerl 
I2009T ). and Si (|Lugarol fl999l ; IZinner et alj l2006h . The de- 
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termination of isotopic ratios from stellar spectra requires 
very high quality data, and isotopic determinations are only 
availa ble for a small n umbe r of elements including Li, C 



(e.g.. ISpite et al.l l200fj) O dSmith fc Lambert! ll990T). M g 
llYong, Lambert, fc Ivans! 120031; iMelendez fc Cohen! \200H ), 
Ti dHughes et alj|2008h . Ba i|Gallagher et alJ^oipfi Tand Eu 



jSneden et al.l l2002). It should be noted however that from 
that list only the carbon isotope ratio is relatively easy to 
obtain. For this reason there is data for the carbon isotope 
ratio for stars in many different evolutionary states, metal- 
licities, and locations (including in external galaxies such as 
the Large and Small Magellanic Cloud). For this reason, the 
13 C/ 12 C ratio is used to study the internal mixing and evolu- 
tion of the observed stars since this ratio changes with time. 
These changes are predicted for low and intermediate-mass 
stars as well as for massive stars that evolve through the 
Wolf-Rayet phase. With the next generation of telescopes, 
the study of isotopes will expand beyond the Milky Way 
and neighbouring Magellanic Clouds toward the brightest 
stars in the outer halo or in neighbouring dwarf spheroidal 
galaxies. This will allow us to gain deeper insights into the 
galactic archaeology beyond our own local neighbourhood. 

However, there are only a small number of the- 
oretical predictions of the evolution of the isotopes 
includ i ng the studies by iTimmes. Wooslev. fc Weaver! 
d 19951). iRomano fc Matteuccil i2003h. IFenner et al.l l|2003L 
IChiappini et alj |2006l ). and lHughes et al.l l|2008l ). In this pa- 
per we present predictions for the time evolution of isotope 
ratios of elemental abundances (from C to Zn) in the Milky 
Way Galaxy. We do this by using the most up-to-date nucle- 
osynthesis yields of core-collapse supernovae and AGB stars 
that are available. In §2, we describe our chemical evolution 
models, and present new and updated supernovae yields. 
These new yield s are then compared to the AGB yields from 
iKarakas! (|2010h . In §3, we show the time/metallicity evo- 
lution of elemental abundances and isotope ratios for the 
solar neighbourhood, bulge, halo, and thick disk. We focus 
on the average evolution of abundances in the Galaxy and 
assume that the ab undances of the carbon-en hanced metal- 
poor stars (CEMP, Bee rs fc Christ lieb 2005]) are explained 
with other effects such as inhomogeneous enrichment, faint 
supernovae, and binary effects. We end with conclusions in 
§4. 



2 THE MODEL 

2.1 Chemical Enrichment Sources 

We include the latest chemical enrichment input into our 
chemical evolution models as follows. The basic equa- 
tions of galactic chemical evo l ution are described in 
Kobayash i. Tsuiimoto fc Nomotol (|2000l . hereafter K00) and 
K06, where the instantaneous recycling approximation is not 
applied, i.e., the contributions from stars of all mass ranges 
are computed as a function of the initial masses and metal- 
licities of stars. 

Stellar winds — The envelope mass and pre-existing 
heavy elements are returned by stellar winds from all dying 
stars. From stars with initial masses of M < 1Mb, it is 
assumed that no new metals are ejected, and that only the 
material outside the He core is returned to the interstellar 



medium (ISM), which contains elements with the abundance 
pattern of the Galaxy at the time when the stars formed. 
The He core mass is set as M re mnant = 0.459 and O.473M0 
for 0.7 and 0.9Mq, respectively. 

Asymptotic Giant Branch (AGB) stars — Stars 
with initial masses between about 0.8 — 8Mq (depending on 
metallicity) pass through the thermally-pulsing AGB phase. 
The He-burning shell is thermally unstable and can drive 
mixing between the nuclearly processed core and envelope. 
This mixing is known as the third dredge-up (TDU), and 
is responsible for enriching the surface in 12 C and other 
products of He-burning, as well as elements heavie r than 
Fe produced by the slow neutron capture process l|Bussd 
lGallind ; lHerwig||2005l ). Importantly, the TDU can result in 
the formation of a C-rich envelope, where the C/O ratio in 
the surface layers exceeds unity. In AGB stars with initial 
masses > AMq , the base of the convective envelope becomes 
hot enough to sustain proton-capture nucleosynthesis (hot 
bottom burning, HBB). HBB can change the surface com- 
position because the entire envelope is exposed to the hot 
burning region a few thousand times per interpulse period. 
The CNO cycles operate to convert the freshly synthesized 
12 C into primary 14 N, and the NeNa and MgAl chains may 
also operate to produce 23 Na and Al. AGB stars with HBB 
have short lifetimes (r < 100 Myr) and are one of the stel- 
lar sites proposed as the p o lluters of glob ular clusters (e.g 
ICottreU fc Da Costal Il98ll ; iRenzinil 120081 ), ev en if quant 



tativ e problems with the models exist (e.g., IFenner et al 



|2004) . Overall a large fraction of light elements such as C 
N and F are produced by AGB stars, while the contribu- 
tion toward heavier elements (from Na to Fe) is negligible 
except perhaps for specific isotopes (e.g., 22 Ne, 25 ' 26 Mg) 
in the context of galactic chemical evolution. AGB stars 
are also an importan t sour ce of elements heavier than Fe 
l|Travaglio et alj|200ll . l2004 ). 

The nucleosynthesis yields of iKarakasl (|2010l ) were cal- 
culated from detailed stellar models, where the structure 
was computed first and the nucleosynthesis calculated af- 
terward using a post-processing algorithm. Yields are in- 
cluded for 77 nuclei including all stable isotopes from H to 
34 S, and for a small group of Fe-peak nuclei. The details of 
this procedure and the codes used to compute the models 
hav e been previou s ly des cribed in some detail, see for exam- 
ple IKarakas et all |2009) and references therein. All mod- 
els were evolved from the zero-age main sequence to near 
the tip of the thermally pulsing AGB. The TDU efficiency 
governs the nucleosynthesis in the lower mass models; this 
was found to vary as a function of H-exhausted core m ass, 
metallicity, and envelope mass (see IKarakas et al.ll2002l . for 
details). For example, in the Z = 0.02 models, no TDU was 
found for M 2Mq . For the intermediate-mass models, the 
TDU was found to be efficient and the occurrence of HBB 
also played a strong role in determining the final yields. The 
occurrence of HBB also strongly depends on the initial mass 
and metallicity, with HBB occurring in lower mass stars with 
a decrease in metallicity (at 3Mq at Z = 10 -4 whereas it 
only starts at ~ 5Mq at Z = 0.02). Furthermore, HBB is 
eventually shut off by the action of mass loss. 

The main uncertainties affecting the nucleosynthesis 
yields of AGB stars involve convection and mass loss. The 
models employ the mixing-length theory of convection with 
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a = 1.75. On the first giant branch, Reimer's mass loss is 
adopt ed with r/B. = 0.4. On the AGB, IVassiliadis fc Wood! 
l)l993l ) mass loss is used for most models, with the exception 
of the intermediate-mass 3 to 6Mq with Z = 10~ 4 mod els, 
where we adopt Reimer's mass loss (see IKarakasl l201(il . for 
details) . 

The m ain differ e nce b etween the iKarakas fc Lattanzid 

\200H ) and IKarakasl ([20101 ) yields is the ch oice of reaction 
rates used in the post-processing algorithm. iKarakasl (|2010T ) 
used an updated set of proton and a-capture rates that 
include some of the latest experimental results for impor- 
tant reactions involve d in the CNO c ycle, NeNa and MgAl 
chains. Furthermore, IKarakasl (|2010h assumed scaled-solar 
initial abundanc es for the Z = 0.008 and Z = 0.004 mod- 
els. In contrast, IKarakas fc Lattanzid (|2007l ') adopted ini- 
tial a bundances for the Large and Small Magellanic Clouds 
from lRussell fc Dopital (Il992l ) which are sub-solar for C, N, 
and O. The updated reaction rates of the 22 Ne(p,7) 23 Na, 
23 Na(p, 7) 24 Mg, and 23 Na(p, a) 20 Ne reactions result in ~ 6 
to 30 times less Na is produced by intermediate-mass mod- 
els with HBB. Note that with t he updated yi e lds th e Na 
overproduction problem found by iFenner et alj (|2004| ) may 
be solved (but not the O depletion r equired stars in glob- 
ular c lusters, although see models bv lVentura fc D'Antonal 
(|2009D ). 

We take the AGB yields for M = 1.0, 1.25, 1.5, 1.75, 1.9 
~ 2.25, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, and 6-5Mpfl 
and Z = 0.0001, 0.004, 0.008, and 0.02 from lKarakasl (|2O10h 7 
We use the yields for Z = 0.02 at Z ^ 0.02. The yields of 
the radioactive isotopes 26 Al and 60 Fe are added to 26 Mg 
and 60 Ni yields, respectively . The masses o f white dwarfs, 
M remnan t, are also taken from lKarakasl ([2010) , and the mate- 
rial outside the white dwarfs (Mi n iti a i — M romnant ) is returned 
to the ISM via stellar winds, which contains the newly pro- 
duced metals (processed metals) and the initial metals that 
existed in the progenitor stars (unprocessed metals). 

The AGB yields (rrn) are defined as the difference be- 
tween the amount of the species in the wind and the ini- 
tial amount in the progenitor star. Therefore, the yields of 
some isotopes (e.g., 15 N) are negative because they are de- 
stroyed during stellar evolution. However, in chemical evo- 
lution models, it is possible that the abundances of such 
elements at time t is lower than the adopted initial abun- 
dances, which causes numerical problems. This is often the 
case for 15 N, and we set the 15 N yield to be if it is negative. 
In chemical evolution models, we define the mass of ejecta 

62 N; 

as the summation of processed metals M j octa = 2J _2 H rrii. 
The rest (M init j a i - M romna nt -M c j ccta ) contains unprocessed 
metals, of which the abundance pattern is not scaled with 
the solar abundance, but is the abundance pattern of the 
galaxy at the time when the stars formed. 

For Z = 0, theoretical models of stars undergo vio- 
lent evolutionary episodes not seen at higher metallicities. 
The ingestion of hydrogen leads to an H flash, followed 
by a "normal" He-shell b urning phase. We take the yields 
and remnant masses from lCampbell fc Lattanzid (|2008T ) for 
M = 0.85, 1.0, 2.0, and 3.0M o . For Na, we assume that the 
yield is reduced by a factor of 10 because the old reaction 



rates were adopted in the calculations. At M > 3Mq, no 
metals are produced. This assumption may not be valid but 
does not affect the average chemical evolution of galaxies. 

Super AGB stars — The fate of stars with initial 
masses between about 8 — IOM0 is uncertain. The upper 
limit of AGB stars, M Ul i, is defined as the minimum mass 
for carbon ignition, and is estimated to be larger at high 
metallicity and also at low metal licity than at Z ~ 10 -4 
(|Gil-Pons et all 120071 ; ISiessI l2007h . At M > M u ,i, stars 
may produce some heavy elements and may explode as so- 
called Type 1.5 supernovae, although no such supernovae 
has yet been observed (K06). Above this ma ss range, stars 
may explode as electron-capture sup ernovae (|Nomotdll984l ; 
iKitaura. Janka. fc Hillebranddlioog 'l . but the metal produc- 
tion (lighter than Fe) is predicted to be very small. In our 
models, we set M Uj i = 4,6.5,6.5,6.5, and 7Mq for Z = 0, 
0.0001, 0.004, 0.008, and 0.02, respectively, and assume that 
no metals are produced from M Uj Jf]to 10M©. The remnant 
mass A:/r Cmnan t is se t as 1.01, 1.12, and 1.15Mq for 7,8, and 
IOMq, respectively. 

Core-collapse supernovae — Although a few groups 
have presented feasible calculations of exploding 10 — 25M© 
stars (|Marek fc Jankall2009l ; iBruenn et al.ll2009H . the explo- 
sion mechanism of core-collapse supernovae (Type II, lb, 
and Ic Supernovae) is still uncertain. However, the ejected 
explosion energy and 56 Ni mass (which decays to 56 Fe) 
can be directly estimated from the observations, i.e., from 
the light curve and spectral fitting of individual super- 
nova. As a result, it is found that many core-collapse su- 
pernovae (M ^ 20M©) have an explosion energy that is 
more than 10 times that of a regular supernova (£51 > 10, 
iNomoto et al.|[2 006j). as well as producing more iron and a 
elements (O, Mg, Si, S, Ca, and Ti). These are called as 
hypernovae (HNe) . The fraction of HNe is uncertain and we 
set e H N = 0-5 at M ^ 20M o . 

iKobavashi et al] (|2006T ) presented the nucleosynthesis 
yields of SNe II and HNe as a function of the progenitor 
mass (M = 13, 15, 18, 20, 25, 30, and 40M Q ) and metallicity 
(Z = 0, 0.001, 0.004, and 0.02). In terms of isotope ratios, 
three models, [M = 18M Q , Z = 0.004] and [M = 25M , 
Z = 0.02] for SNe II and HNe, showed a relatively large 
production of 13 C and N. This was caused by the convective 
mixing of hydrogen into the He-burning layer. The cause 
of this mixing is uncertain and it does not always manifest 
itself in the stellar models. For this reason, we re-calculated 
the progenitor star models and explosive nucleosynthesis. 
The results of these new calculations without such mixing 
are presented in Tableland shown in Figs. IHlOl While the 
updated yields of 13 C and N are reduced compared to the 
results of K06, the yields of other major isotopes are not 
significantly different. 

As in K06, the yield tables provide the amount of pro- 
cessed metals (m;) in the ejecta (in Mq). The mass of the 

ejecta is given as M ojccta = X^=?H m * = ^flnai - Af cu t- 
Stellar winds reduce the stellar mass to M nna i at the on- 
set of the supernova explosion, with the central mass M cu t 
falling onto the remnant. The stellar winds with the mass of 
A/initi a i — Affl na i contains unprocessed metals, which are not 
included in the tables but are included in the chemical evo- 



1 2.1M model is provided for Z = 0.004 and 0.008. 

2 6.5Mq model is available only for Z = 0.02. 



3 The remnant mass at M Ui i is extrapolated for Z < 0.02. 
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lution models. The contribution from the stellar winds are 
added (Eq.8 in K00) to the supernova ejecta (Eq.9 in K00). 
The abundance pattern of the stellar winds is not scaled 
with the solar abundance but is the abundance pattern of 
the galaxy at the time when the stars formed. Note that 
newly synthesized 4 He in the winds is not included and thus 
we will leave the evolution of helium to a future study. 

The upper limit of core-collapse supernovae, M u ,2, is not 
well known owing to uncertainties in the physics of blackhole 
formation. We set M Ui 2 = 50Mq, which is constrained from 
the [a/Fe] plateau at [Fe/H] < - 1 (K06). Because the 
envelope mass that contains a elements is larger for massive 
stars, a larger M Ut 2 results in higher [a/Fe] ratios. For M > 
M n< 2, we assume that no metals are produced because the 
central part of the progenitor is likely to collapse to form a 
blackhole. 

Rotating massive stars — The rotation of stars 
induces the mixing of C into the H-burning shell, 
which produces a larg e amount of primary nitrogen 
jMevnet fe Maederl 12003 : Iffirschil I2007T ). Rotation also af- 
fects the CO core mass, and a few groups are calculating 
the stellar evolution of rotating stars from the main se- 
quence through to the final explosive supernova stage, al- 
though the yields are not yet available. The rotational ve- 
locity of these stars is typically a free parameter. To show 
the effect of rotation, in one of our models we incl ude the 
yields of 3 ' 4 He, 12 ' 13 C, 14 N, 16 - 17 ' ls O, and 22 Ne from lHirschil 
(|2007D for Z = with [20M©, 600 km s~l], [4OM , 700 km 
s"l], and [85M , 800 km s"l], with M u>2 = 100M Q . 

Pair-Instability supernovae — Stars with WOMq < 
M < 300Mq encounter the electron-positron pair insta- 
bility and do not reach the temperature of iron photodisin- 
tegration. Pair-instability supernovae (PISNe) are predicted 
to produce a large amount of metals such as S and Fe. This 
contribution is not included in our chemical evolution mod- 
els because the number of such stars are expected to be 
very small and because no signature of PISNe has been de - 
tected in metal-poor stars dUmeda fe Nomot3 20021 . l2005h . 
Recently, iKobavashi. Tominaga. fe Nomotol ( 2010T ) showed 
that the observed elemental abundance pattern of very 
metal-poor Damped Lyman a systems is inconsistent with 
the yields of PISNe, and is instead consistent with the yields 
of faint core-collapse supernovae. 

Type la Supernovae (SNe la) — The progenitors 
of the majority of SNe la are most likely Chandrasekhar 
(Ch) mass white dwarfs (WDs). For the evolution of ac- 
creting C+O WDs toward the Ch mass, two scenarios have 
been proposed. One is the double-degenerate scenario, i.e., 
the merging of double C+O WDs with a combined mass 
surpassing the Ch mass limit. However, it has been theo- 
retically suggested that t his leads to accretion-ind uced col- 
lapse rather than SNe la l|Nomoto fe K ondo 199l]), and the 
lifetimes are too short to repr oduce the chemical evo lution 
in the so lar neighbourhood ( Kobavas hTet al.l 1 19981 . here- 
after K98: lKobavashi fe Nomotol2009l . hereafter KN09). The 
other is the single-degenerate (SD) scenario, i.e., the WD 
mass grows by accretion of hydrogen-rich matter via mass 
transfer from a binary companion. The mass accretion rate 
is limited to trigger carbon deflagration (|Nomotdll983 '), but 
the allowed parameter space of binary systems can be sig- 
nificantly increased by the WD wind effect if the metallicity 
is higher than [Fe/H] ~ -1 (K98). 



30 



20 



X 



10 



1 1 1 1 1 1 1 II 
- z = o 


i i i i i i i 










13 C 




- 


14 N 






15 N 






16Q 






18Q 




_ 


19JP 






25 Mg 






a6 Mg 


/ / 




56 Fe 


• — ~~ / 














10 



20 

M [M ] 



30 



40 



Figure 1. Mass fraction elements in the ejecta for Z = as a 
function of progenitor ma ss, normalized to the solar abundances 
llAnders fe G revessc 198^). The solid and dashed lines are for SNe 
II and HNe, respectively. 
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In our models, based on the SD scenario, the lifetime 
distribution function of SNe la is calculated with Eq.[2] 
in KN09, taking into account the metallicity dependence 
of the WD winds (K98) and the mass-stripping effect on 
the binary companion stars (KN09). There are two kinds 
of progenitor systems. One is the main-sequence+WD sys- 
tem with timescales of ~ 0.1 — 1 Gyr, which are dominant 
in star- forming galaxies (the so-called prompt population). 
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Figure 4. The same as FigQ]but for Z = 0.008. 



The other is the red-giants+WD system with lifetimes of 
~ 1 — 20 Gyr, which are dominant in early-type galaxies. 
Although the metallicity effect of SNe la has not yet been 
confirmed by supernova surveys, it is required to account for 
the presence of a young population of SNe la, which in turn 
are required by chemical evolution of the Milky Way Galaxy 
(KN09). Note that the observed elemental abundance pat- 
tern (e.g., the low [Mn/Fe] ratios) in dwarf spheroidal galax- 
ies are more consistent with the enrichment from low-mass 
SNe II than that from the SN la enrichment (K06). 

For SNe la we take the nucleosynthesis yields from 



iNomoto et al.l (|l997h . The metallicity dependence of the 
progenitors is not included but is not expected to be very 
large (H. Umeda, K. Nomoto, et al., private communica- 
tion). Note that Ni is overproduced at [Fe/H] > — 1 com- 
pared to the observations in the solar neighbourhood, but 
this can be solved by tuning the propagation speed of the 
burning front and th e central density of the white dwarf 
(|lwamoto et al.lll999h . 

2.2 Yields 

Figures [T][S] and Figures IfJl llOl show the yields and the isotope 
ratios as a function of t he mass and metallicity, n ormalized 
to the solar abundance (|Anders fe Greve sse 19891 hereafter 
AG89). The yields are also normalized to the mass of the 
ejecta M c j ec ta, which are much smaller for AGB stars than 
for supernovae. However, the AGB yields of C, N, and F are 
comparable to those of supernovae. 

CNO cycling at T > 2 x 10 7 K results in the pro- 
duction of 13 C, 14 N, and 17 0. Core and shell He-burning 
results in the synthesis of 12 C, ie O, and 19 F at T > 
1.5 x 10 8 K. The heavy isotope of oxygen, ls O, is de- 
stroyed by proton captures in stellar interiors by CNO cy- 
cling, and produced and then destroyed again by a-captures 
via 14 N(q|7 ) 18 F( / 9 + ) 18 followed by 18 0(«, 7 ) 22 Ne (e.g., 
lArnettl [l996l ). Further secondary He-burning reactions can 
synthesize 25 ' 26 Mg via 22 Ne(a, n) 25 Mg and 22 Ne(a, 7) 26 Mg, 
which follow after ls O(a, 7 ) 22 Ne. The majority of 18 and 
elements heavier than fluorine up to the Fe peak are pro- 
duced by hydrostatic burning (T > 7 x 10 8 K) and explosive 
nucleosynthesis (T >2x 10 9 K) in massive stars (M >8M ). 

The AGB yields show that 12 C and 19 F are signifi- 
cantly produced in low-mass stars (1 — 4Mq) whereas 13 C, 



14 N, 25 Mg, and ^Mg are produced in intermediate -mass 
(4 - 7M Q ) stars. The yields of 19 F in iKarakasj (|2010| ) have 

owing 



21) » 



increased compared to iKarakas fc Lattanziol 1 20071 ) 
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Figure 6. Mass ratios of isotopes in the ejecta for Z = as 
a function of progenitor mass, normalized by the solar ratios 
llAnders fc G rcvcssc 1989) in the logarithmic scale. The solid and 
dashed lines are for SNe II and HNe, respectively. 



to a slower 19 F(a,p) 22 Ne reaction rate. The production of 
F and the neutron-rich Mg isotopes in AGB stars are highly 
mass depen dent. F production peaks at ~ 3M Q at solar 
metallicity ()Lugaro et al.l[2004h ; in higher mass models F is 
destroyed by a-captures by the higher temperatures reached 
during He-burning. Likewise, 25 Mg and 26 Mg are only pro- 
duced by He-burning when the temperature exceeds about 
~ 300 x 10 6 K and these conditions are only reached by 
intermediate-mass AGB stars karakas et al. l B). "O is 
produced by the entire AGB mass range, with relatively 
more from low-mass AGB stars. The yields of 15 N, ls O, and 
18 O are negative in most cases indicating that these isotopes 
are destroyed by AGB nucleosynthesis with the largest de- 
struction taking place by HBB in intermediate-mass AGB 
stars. 

Core-collapse supernovae are the main producers of the 
major isotopes with more minor isotopes produced at higher 
metallicity. This is because minor isotopes are synthesized as 
secondary elements from the seed of the major isotopes. As a 
result, the ratios between the major and minor isotopes are 
larger for supernovae than for AGB stars, except for 12 C/ 13 C 
at 1-4Mq (and 14 N/ 15 N). The ratios are, in general, larger 
for low-metallicity supernovae, and approach the solar ratios 
([ 1 X/ 2 X]=0) with an increase in the metallicity. 



2.3 Initial Mass Function 

We adopt the recent obs ervational esti mate of the initial 
mass function (IMF) from lKro upa (200§|), which is a power- 
law mass spectrum 4>(m) oc m~ x with three slopes at dif- 
ferent mass ranges: x = 1.3 for 0.5Mq m ^ 50Mq, 
x = 0.3 for O.O8Af sC m < 0.05M©, and x = -0.7 for 
0.01M Q < m < 0.08M Q . At the high-mass end, the s l ope o f 
the Kroupa IMF is almost the same as the lSalpeten (|l955t ) 



Figure 7. The same as FigE]but for Z = 0.001. 
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Figure 8. The same as Fig[6]but for Z = 0.004. 



IMF (x = 1.35), and is flatter than the lMiller fc Scald (|l979h 
IMF. Metal enrichment is only obtained from stars with m > 
0.5Mq, hence chemical evolution predictions obtained with 
the Kroupa IMF are not significantly different from K06's 
results with the Salpeter IMF at 0.07M Q < m < 50M o . In 
Figure[TT]we show that the Kroupa (short-dashed lines) and 
Salpeter (solid lines) IMFs can provide almost the same age- 
metallicity relation, [0/Fe]-[Fe/H] relation, and metallicity 

distribution function (MDF). 

The mass distribution of stars from the lChabriej ([2003) 
IMF is peaked at ~ 3Mq, leading to fewer stars with 
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m < 0.5Mq, and a steeper slope at the high- mass end than 
the Kroupa and Salpeter IMFs. Although the Chabrier IMF 
is adopted for many cosmological simulations, the metal pro- 
duction that results is too large to meet the observational 
constraints of the s olar neighbourhood. In particular the 
present SN la rate dMannucci et al]|2005l ) is too high and 
the [a/Fe] is too high as shown in Figure [TT] (long-dashed 
lines). These issues cannot be solved by changing SN la pa- 
rameters. 



Table 1. Parameters of chemical evolution models: Infall, star 
formation, and outflow timescales, and the galactic wind epoch 
in Gyr. 
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2.4 Star Formation Histories 

We use similar models as K06 for the star formation histories 
of the solar neighbourhood, bulge, halo, and thick disk, but 
with the Kroupa IMF. The galactic chemical evolution is 
calculated with the basic equations described in K00 and 
K06. In one-zone chemical evolution models, the gas fraction 
and the metallicity of the system evolve as a function of time 
by star formation, as well as infall and outflow of matter. The 
star formation rate (SFR) is assumed to be proportional to 
the gas fraction (— / g ). The infall of primordial gas from the 
outside of the component is given by the rate oc texp[— 
for the solar neighbourhood, and ^exp( — for the other 
components. For the halo, outflow is included such that the 
rate is proportional to the star formation rate ( — / g ). For 
the bulge and thick disk, star formation is assumed to be 
truncated by galactic winds at a given epoch (t = t w ). 

The MDF is one of the most stringent constraints on 
the star formation history. The parameters that determine 
the SFRs are chosen to meet the observed MDF of each 
component as summarized in Table [1] The MDFs, the re- 
sultant SFR histories, and the age-metallicity relations are 
shown in Figure [T2j The uncertainties of the observations 
were discussed in K06 and will not be repeated in this pa- 
pe r. The observed MDF in the Galactic bulge was updated 
by IZoccali et al.l (120081 ) . where the peak met allicity is sig- 
nifican tly higher than the previous MDF by IZoccali et al.l 
(2003). Except for the bulge model, compared with the K06 
results, the small differences originate from the choice of the 
IMF and not from the updated yields presented in this pa- 
per. 

In the solar neighbourhood model (solid lines), star for- 
mation takes place over 13 Gyr. The MDF shows a narrow 
distribution peaked around [Fe/H] ~ —0.2, which is consis- 
tent with the observations. Note that introducing infall sig- 
nificantly redu ces the number of metal-poor G-dwarf stars 
(|Tinsley|| 19801 ) , and thus there is no G-dwarf problem in the 
predicted MDF according to our models. 

For the bulge (long-dashed lines), we use the 
infall+wind model (model B in K06) with a short star- 
formation timescale. The infall is required to explain the lack 
of metal-poor stars and the wind is adopted to reproduce the 
sharp-cut of the observed MDF at the metal-rich end. It is 
possible that star formation continues at present time, form- 
ing super metal-rich stars in the Galactic bulge. In this case 
the SFR should be low to meet the MDF. In our bulge model 
the duration of star formation is set to be 3 Gyr, which re- 
sults in a peak metallicity of [Fe/H] ~ +0.3. A much higher 
efficiency of chemical enrichment, e.g., a flatter IMF is not 
required, unless the duration is much shorter than 3 Gyr. 
Note that the 3 Gyr duration is consistent with chemody- 
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namical simulations of Milky Way-type galaxies from CDM 
initial conditions (KN11). 

For the thick disk (dot-dashed lines), we use the 
infall+wind model (model C in K06), which also gives a 
;ood agreement wit h the observed age-metallicity relation 
Bensbv et al]|2004) . In the thick disk model, the formation 
timescale is as short as ~ 3 Gyr, and the star formation 
efficiency is larger than that for the solar neighbourhood 
but smaller than for the bulge. The assumption that the 
timescale of star form ation is shorter in the bulge a nd thick 
disk was suggested by iMatteucci fc Brocatol (ll990T ). but by 
itself does not completely explain the observed abundance 
patterns. A more intense star formation is also required as 
well as a shorter star formation timescale. 

For the halo (short-dashed lines), we use an outflow 
model without infall, wh ich results in a low me an metallicity 

of [Fe/H] 1.6 (e.g.. IChiba fc Yoshiil Il998h . In the halo 

model, the star formation efficiency is much lower than for 
the other components, and the outflow causes an effective 
metal loss. 

The me tallicity dependent m a in-seq uence lifetimes are 
taken from iKodama fc Arimotol (|l997h for 0.6 - 8OM , 
which are calculated with the st ellar evolution code de- 
scribed in llwamoto fc Said dl999h. These are in excellent 
agreement with the lifetimes in iKarakas! (|2010h for low- and 
intermediate-mass stars. 



3 RESULTS 

Figure [T3] shows the evolution of element abundance ratios 
[X/Fe] against [Fe/H] for the solar neighbourhood, using the 
models with SNe II, HNe, and SNe la only (dashed lines), 
with AGB stars (solid lines), and with rotating massive stars 
at Z — (dotted lines) . The difference among models is seen 
only for C, N, and F. For the elements heavier than Na, the 
small difference among the models with and without rotat- 
ing massive stars is caused by the difference in the IMF 
(M u ,2 = 12OM0 is adopted in the case with rotating mas- 
sive stars instead of M Ui 2 = 50Mq). Observational data are 
taken from several sources, which were selected to minimize 
systematic errors as discussed in KNU. In that paper, we 
also discussed the different results for the lines used in the 
observational data analysis, which will not be repeated in 
this paper. For C and N, only data from unevolved stars are 
plotted. Both in the models and the observational data, the 
solar abundances from AG89 are adopted. The results are 
summarized as follows. 

• a elements — In the early stages of galaxy formation 
only SNe II/HNe contribute and the [a/Fe] ratio quickly 
reaches a plateau ([a/Fe] ~ 0.5). Around [Fe/H] ~ —1 SNe 
la start to occur, which produce more iron than a elements. 
This delayed enrichment of SNe la causes the decrease in 
[a/Fe] with increasing [Fe/H]. The [Fe/H] where the [a/Fe] 
starts to decrease depends on the adopted SN la progeni- 
tor model, and is determined not by the lifetime but by the 
metallicity dependence of SN la progenitors (KN09). As a 
result, this trend is in excellent agreement with the obser- 
vations (dots) for O, Mg, Si, S, and Ca. Ne and Ar also 
show a similar trend. Ti is underabundant overall, but this 
problem can be solved w ith 2D nucleosynthesis calculations 
dMaeda fc Nomotoll2003h . as well as for Sc and V. AGB stars 



do not make any difference to these trends. Although AGB 
stars produce significant amounts of the Mg isotopes, the 
inclusion of these do not affect the [Mg/Fe]-[Fe/H] relation. 

• Odd-Z elements — Na, Al, and Cu show a decreasing 
trend toward lower metallicity, which is well reproduced by 
the strong dependence of these elements on the metallicity of 
the progenitor stars (see Fig. 5 in K06). In contrast, Na and 
Al show a decreasing trend toward higher metallicity owing 
to the contribution from SNe la, which is shallower than the 
trend for the a elements. Such a decrease is not seen for 
Cu since Cu is also produced by SNe la. With the updated 
AGB yields (§2.2), [Na/Fe] is consistent with the observa- 
tions, and the Na overproduction problem by AGB stars is 
not seen. Note that AGB stars may produce some Cu but 
no yields are available for a large range of masses and metal- 
licities. K, Sc, and V are underabundant overall, a problem 
which has not been discussed in detail in previous studies. 
The i/-process can increase the production of these elements 
(Izutani, Umeda, & Yoshida 2010, private communication), 
although the yields are not yet available. [(P, Cl)/Fe] are 
also negative overall in our predictions. There is a metal- 
licity dependence of P, CI, K, and Sc yields at Z > 0.001 
for SNe II/HNe, which causes a weak decrease from [Fe/H] 
~ — 1 to ~ —3. The V yields do not depend very much on 
metallicity. 

• Iron-peak elements — [(Cr, Mn, Co, Ni, Zn)/Fe] are 
consistent with the observed mean values at —2.5 < [Fe/H] 
<^ — 1. Note that Cr II observations are plotted, be- 
cause this line is not strongly affected by NLTE effect. For 
Mn, the NLTE effect sh ould not be so large as indicated 
by M n II observations (|johnsonl |2002| ; iMashonkina et all 
2010), although a strong NLTE effect is reported by 
Bergemann fc Gehrenl (|2008T ). At [Fe/H] < - 2.5, observa- 



tional data show an increasing trend of [(Co, Zn) /Fe] toward 
lower metallicity, which will not discussed here since inhomo- 
geneous chemical enrichment is becoming increasingly im- 
portant. The [(Co, Zn)/Fe] trend can be explained by HNe 
under the assumption that the observed stars wer e enriched 
by only a single supernova (|Tominaga et al.ll2007l ). 

• Manganese — Mn is a characteristic element of SN 
la enrichment and is produced more by SNe la than 
SNe II/HNe relative to iron. From [Fe/H] ~ -1, [Mn/Fe] 
shows an increasing trend toward higher metallicities, 
which is caused by the del a yed e nrichment of SNe la. 
iFeltzing. Fohlman. fc Bensbvl |2007T) showed a steep slope 
at [Fe/H] > 0, which could b e generated by the m etallicity 
dependence of SN la yields. ICescutti et al.l (|200ST ) demon- 
strated this by applying an artificial metallicity dependence 
to the theoretical calculations. In principle, Mn is an odd-Z 
element and the Mn yields depend on the metallicity both 
of SNe II and SNe la. The metallicity dependence for SNe II 
are included in our yields. For SNe la, a strong metallicity 
dependence is not expected (H. Umeda, K. Nomoto, et al., 
private communication) . 

• Zinc — Zn is one of the most important elements for 
supernova physics. [Zn/Fe] is about ~ for a wide range 
of metallicities, which can only be generated by a large 
fraction of HNe (50% of M ^ 20M©). In detail, there is 
a small oscillating trend; [Zn/Fe] is at [Fe/H] ~ 0, this 
increases to 0.2 at [Fe/H] ~ —0.5, decreases to be again at 
[Fe/H] ~ —2, then increases toward lower metallicity. This 
is characteristic of our SN la model (KN09) and is consis- 
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Figure 11. Star formation histories (panel a), age-metallicity relations (panel b), [0/Fe]-[Fe/H] relations (panel c), and metallicity 
distribution functions (panel d) for the solar neighbourhood, with the Kroupa (2008) IMF (solid lines), the Salpeter ( 1955) IMF ( short - 
dashed lines), and the Ch abrier (2003) IMF (long- dashed lines). The observational da ta sources are: an error e stimate, iMatteuccil l|l997l l 
in panel (a); filled circ les. | E dvardsspn et ah! lll993l ) in pan els (b) and (d); open circles, I Wvse &; G ilmorel Jl995h in p a nel (d ). In the panel 
(c), large filled circles. ICavrel et al.l (j2004T ): filled squares. IcTratton et all <|2003h : filled triangles. iBensbv fc Feltzingl j2006t) . 



tent with the observations l|Saito et al.ll2009T ). Theoretically, 
Zn production depends on many parameters; 64 Zn is synthe- 
sized in the deepest region of HNe, while neutron-rich iso- 
topes of zinc 66_70 Zn are produced by neutron-capture pro- 
cesses, which are larger for higher metallicity massive SNe II. 
Since the observed [Zn/Fe ] ratios show an increasing trend 
towa r d lower metallicit y (|Primas et al.l l200d ; iNissen et al.l 
120071 ; ISaito et al.1 120091 ). the HN fraction may have been 



larger in the earliest stages of galaxy formation. At higher 
metallicities, the HN fraction may be as small as 1% (KN11). 

• Carbon — Although the ejected mass of C is similar 
for low-mass AGB (1 — 4Mq) and massive (> IOMq) stars, 
the [C/Fe] ratio is enhanced efficiently by low-mass stars 
because these stars produce no Fe (Figs. IU5[I . When we in- 
clude AGB yields (solid lines), [C/Fe] increases from [Fe/H] 
~ —1.5, which corresponds to the lifetime of ~ 4Mq stars 
(~ 0.1 Gyr). At [Fe/H] 1, [C/Fe] reaches 0.13, which is 
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Figure 12. Star formation histories (panel a), age-metallicity relations (panel b), and metallicity distribution functions (panel c) for the 
solar neighbourhood (solid lin es), halo (short-dashed li nes) , bulg e (long-dashed lines) , a nd thick disk (dot-dash ed lines). The observational 
data sources are: filled c ircles, [Ed vardsso n et alj l|l993l l: crosses, IChiba fc Yoshiil l|l99Sl ): filled triangles. IZoccali et alj (|2008h : open circles, 
IWvse fc Gilmord il995l) . 



0.32 dex larger than the case without AGB yields (dashed 
li nes). This is roughly consistent with previous m odels such 
as lPrantzos. Vangioni-Flam. fc Chauveaul (| 19941 ). At [Fe/H] 
> — 1, [C/Fe] shows a decrease due to SNe la. Because 
of the long lifetimes of AGB stars, no difference is seen at 
[Fe/H] < — 1.5, which is consistent with the observed be- 
havior of s-process elements (|Travaglio et al.ll2004l ) . If we in- 
clude the yields of rotating massive stars (dotted lines), the 
[C/Fe] ratio becomes as large as ~ 0.5 at [Fe/H] ~ —2.5. 
A significant fraction of metal-poor stars show carbon en- 



richment (CEMP stars), with several scenarios proposed to 
explain t he observed abundances including a single super- 
nova (e.g lllmeda fc Nomotcll2002l) and AGB s t ars in binary 
syste ms l|Suda et alj|2004 iLugaro et all 120081 ; llzzard et all 
120091 ). Such local peculiar effects are not included in our 
models. We should also note that AGB stars can contribute 
at metallicities below [Fe/H] < — 1.5 when an inhomoge- 
neous chemical enrichment is taken into account (KN11). 

• Nitrogen — Different from C, N is produced mainly by 
intermediate-mass AGB stars (4 — 7Mq , independent of the 
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Figure 13. Evolution of elemental abundance ratios [X/Fe] against [Fe/H] for the solar neighbourhood with SNe II, HNe, and SNe la 
only (dashed lines), with AGB stars (solid lines), and with rotating massive stars at Z = (dotted lines). The dots are observational 
data (see KN11 for the references). For C and N, only unevolved stars are plotted. 
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Figure 14. Evolution of elemental abundance ratios [X/Fe] against [Fe/H] for the solar neighbourhood (solid lines), halo (short-dashed 
lines), bulge (long-dashed lines), and thick disk (dot-dashed lines) with AGB yields. The black and yellow dots are observational data 
for thin and thick disk stars (see KN11 for the references). 
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integration on the IMF). Therefore, the contribution from 
AGB stars (solid lines) is seen at [Fe/H] ~ -2.5. At [Fe/H] 
~ — 1, [N/Fe] reaches 0.15, which is 0.8 dex larger than the 
case without AGB yields (dashed lines). At [Fe/H] > — 1, 
[N/Fe] shows a shallow decrease due to SNe la. No differ- 
ence is seen at [Fe/H] < — 2.5 with and without the AGB 
yields, while [N/Fe] can be as large as ~ 0.5 with rotating 
massive stars (dotted lines). Ichiappini" et all (|2006l ) showed 
that the contribution from rotating massive stars is required 
to solve the primary N problem. As noted above, however, 
AGB stars can also contribute to N production even at 
[Fe/H] < — 2.5 when taking inhomogeneous chemical en- 
richment into account (KN11). From the difference between 
C and N, it is possible to distinguish the contribution from 
low- and intermediate-mass AGB stars as the enri c hment 
source of the observed metal-poor stars. iPols et al.l (|2009h 
showed that the IMF with the G aussian distribution peaked 
at ~ lOM© (|Komiva et al.ll2007l ) is rejected when trying to 
match the fraction between C-rich and N-rich stars, based 
on bi nary population synthesis models (see also llzzard et al.l 
l2009h . 

• Fluorine — F is one of the most interesting elements, 
although F abundances are estimated from only one infrared 
line from stellar spectra. AGB stars and massive stars have 
both been suggested to produc e F but production has only 
been confirmed for AG B stars jjorissen. Smith, fc Lambert! 
ll992l ; lAbia et al.ll201ol ). The AGB mass range that produces 
F is similar for C and is 2 — 4Mq. Thus the difference is 
seen only at [Fe/H] > - 1.5. At [Fe/H] ~ -1, [F/Fe] reaches 
0.22 in the model with the AGB yields (solid lines). This 
is 0.56 dex larger than the case without the AGB yields 
(dashed lines) and much closer to the observational data 
|Cunha et alj|2003h . Note that the F yields from AGB stars 
were increased with the new reaction rates (§2). Different 
from C, F is not significantly produced by SNe II/HNe ac- 
cording to our yields, and thus [F/Fe] rapidly decreases from 
[Fe/H] ~ —1 to ~ —3. Therefore, the F abundance is a good 
clock to distinguish the contribution from low-mass AGB 
stars and supernovae. We should note, however, that the 
F yields from supernovae may be increased by a factor of 
~ 1000 by the ^-process (Izutani , Umeda, fc Yoshida 2010 , 
private communication; see also (Wooslev fc Weaverl Il995h . 
The effect of rotating massive stars is uncertain since F 
yields are not available in the literature. 

These results are consistent with the mo dels that adopt 
the s ame supernova and AGB yields in iRomano et al.l 
l|2010l. hereafter R10). There are some differences, for ex- 
ample, for elements heavier than N, our results fall be- 
tween Models 4 and 5 in R10 because we adopt an hy- 
pernova efficient of ehn = 0.5. The C and N abundances 
are predicted to be higher in Models 4 and 5 than our re- 
sults because the low and inte r media te-mass star yields of 
Ivan den Hoek fc Groenewegenl (1 19971 ) were used. In Model 
15 of R10, the same AGB yields are adopted as in our mod- 
els, but the C and N abundances are still higher. This may 
be due to the addition of the Geneva pre-supernova yields 
along with the supernova yields from K06. We note that ro- 
tation changes the pre-supernova structure and thus should 
also change the nucleosynthesis during the explosion. Hence 
to be fully self consistent, supernovae yields computed from 
a rotating pre-supernova structure should be included but 
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Figure 15. Evolution of the [C/O] ratio against [O/H] for the 
solar neighbourhood (solid lines), halo (short-dashed lines), bulge 
(long-dashed lines), and thick disk (dot-dashed lines) with AGB 
yields. The observat i onal d ata sources are: filled circles, unmixed 
stars in ISpite et al.l 1I2OO6I ) ; filled triangles and open triangles, 
iBensbv fc Feltzinsd 1I2OO6I) for thin a nd thick disk stars, respec- 
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Figure 16. The same as Fig. [TS] but for the [F/O] ratio. 
The o bservational data sources are: open circles, ICunha et al.l 
(2003) for the solar neigh bourhood stars; filled circles, 
ICunha. Smith, fc Gibson! teOOcj) for bulge stars. 



14 Chiaki Kobayashi, Amanda I. Karakas, and Hideyuki Umeda 



these are not yet available for a wide range of masses and 
metallicities. Other model differences include the adopted 
IMF, star formation rates, and the SN la model, but these 
do not significantly affect the average evolution of elemental 
abundance ratios. However, we do observe small differences 
in e.g., the [a/Fe]-[Fe/H] relations as a result of the the 
different adopted SN la model. In RIO's models, the evolu- 
tionary change in [O/Fe] around [Fe/H] ~ — 1 is not as sharp 
as in our models. In addition, the evolutionary track around 
[Fe/H] ~ —0.8 in the RIO models is not smooth, which may 
result in a large number of stars at this point in the [O/Fe]- 
[Fe/H] diagram. Large and homogeneous observational data 
sets may help to put constraints on the modelling of SNe la. 

Enrichment sources produce different elements on dif- 
ferent timescales, and thus the time evolution of the elements 
varies as a function of location in a galaxy, depending on the 
star formation history. In Figure [14] we show the evolution 
of elemental abundance ratios [X/Fe] against [Fe/H] for the 
solar neighbourhood (solid lines), halo (short-dashed lines), 
bulge (long-dashed lines), and thick disk (dot-dashed lines), 
where the contributions from SNe II, HNe, SNe la, and AGB 
stars are included. 

Bulge and thick-disk — If the star formation 
timescale is shorter than in the solar neighbourhood (solid 
lines) as in our bulge (long-dashed lines) and thick-disk 
(dot-dashed lines) models, the contribution from stars of 
a given lifetime appear at a higher metallicity than in the 
solar neighbourhood. Intermediate-mass AGB stars, low- 
mass AGB stars, and SNe la start to contribute at [Fe/H] 
~ —2.5, —1.5, and —1, respectively in the solar neighbour- 
hood, but at a higher [Fe/H] in the bulge and thick disk mod- 
els. At [Fe/H] > - 1, [a/Fe] is higher and [Mn/Fe] is lower 
than in the solar neighbourhood because the SN la contribu- 
tion is smaller in the bulge and thick disk. Simultaneously, 
the [(C, N, F)/Fe] ratios peak at higher metallicities; [(C, 
N, F)/Fe] is lower at [Fe/H] < — 1, and is slightly higher at 
[Fe/H] > — 0.5 than in the solar neighbourhood. The abun- 
dance ratios of [(Na, Al, Cu, Zn)/Fe] and [(P, CI, K, Sc)/Fe] 
are predicted to be higher because of the higher metallic- 
ity. Indeed, the metallicity reaches values high enough to 
produce these elements before the majority of SNe la oc- 
cur. The yellow dots are the observational data of the thick 
disk stars, which show higher [a/Fe] and [(Al, Cu)/Fe] ratios 
than the thin disk stars (black dots) . These data are roughly 
consistent with our model predictions. Note that the tran- 
sition metallicity where [a/Fe] starts to decrease and the 
peak metallicity of [(C, N, F)/Fe] and [(Na, Al, Cu, Zn)/Fe] 
depends on the star formation and infall timescales of the 
system. For the bulge model, [O/Fe] starts to decrease at 
[Fe/H] ~ —0.5. The transition met allicity can be increa sed 
with a much shorter timescale as in iBallero et all (|2007h . 

The high [a/Fe] in the bulge and thick disk can also be 
reproduced by changing the IMF, namely, adopting a flatter 
IMF (see Fig. 32 of K06). In this paper we do not need to 
change the IMF between the solar neighbourhood and the 
bulge (§2.4). This is beca use the observationa l constraints of 
the MDF was updated bv lZoccali et al.l (|2008l ') and the bulge 
stars are more metal-rich than in the solar neighbourhood 
(Fig[T2]). With a flatter IMF, the predicted [(Zn,Co)/Fe] ra- 
tios becomes much larger, a result that can be tested with 
future observations. 



The difference in the AGB contribution is clearly seen in 
[(C, N, F)/0] ratios. Figure [Pol shows the evolution of [C/O] 
against [O/H], where [C/O] shows a rapid increase at [O/H] 
~ — 1 to ~ —0.5 toward higher metallicities. At [O/H] > — 1, 
the [C/O] ratio is highest in the halo, followed by the solar 
neighbourhood, thick disk, and bulge. This is qualitatively 
consistent with the observations (dots) in the thin disk, thick 
disk, and bulge, but the very high [C/O] ratios that are ob- 
served in som e stars at [O/H] ~ cannot be reproduced 



by our models. ICescutti et alj (|2009h showed that the [C/O ] 
ra tios could be enhanced a t [O/H] > if iMaederl (I1992T ) 
or iMevnet fc Maed"erl (I2002T ) yields are included. Similarly, 
Figure \W\ shows the evolution of [F/O] against [O/H], where 
[F/O] shows a rapid increase at [O/H] ~ —1.5 to ~ —1 to- 
ward higher metallicities. The present [F/O] ratio is slightly 
lower than the observations, per haps indicating the need for 
other sources of F in the galaxy. iRenda et alj (|2004l ) showed 
that [F/O] could be enhanced at [O/H] > - 0.2 if the yields 
of Wolf-Rayet stars ar e inclu ded. 

iMcWiiliam et all (120081 ) showed that the [O/ Mgj ra- 
tio co uld be changed in their bulge models if the IMaederl 
l|l992D 's yields are included. This ratio cannot be changed 
by the addition of AGB yields. With our yields, [O/Mg] is 
almost constant independent of the metallicity as shown in 
Fig. 9 of K06. Our model predictions are consistent with the 
observations of thin disk stars, but not with those of thick 
disk and bulge stars, although the scatter and uncertainties 
are quite large for the bulge stars. 

Halo — If the chemical enrichment timescale is longer 
than in the solar neighbourhood as in our halo model (short- 
dashed lines) , the contribution from low-mass AGB stars be- 
come significant, and thus the [(C, F)/Fe] ratios are higher 
at all metallicities than in the solar neighbourhood. For 
[N/Fe], the difference is seen only at low metallicity ([Fe/H] 
< — 1) since the mass range of N production is more mas- 
sive than that of C and F. The [a/Fe] and [Mn/Fe] relations 
are the same as in the solar neighbourhood. However, the 
[(Na, Al, Cu, Zn)/Fe] and [(P, CI, K, Sc)/Fe] ratios are lower 
because of the overall lower metallicity. Detailed elemental 
abundances are not available for the Galactic halo, but high 
[C /Fe] is seen for a significant fracti on of stars in the Galac- 
tic outer halo in the SEGUB^ data l|Beersll201fj| '). 

Figure [17] shows the evolution of the isotope ratios 
against [Fe/H] for the solar neighbourhood (solid lines). 
In general, core-collapse supernovae are the main produc- 
ers of the major isotopes with more minor isotopes syn- 
thesized at higher metallicity (see Figs. |6"ITT7)) . For this rea- 
son, the evolution of 12 C/ 13 C, 16 0/ 17 ' ls O, 20 Ne/ 21 > 22 Ne, 
24 Mg/ 25 ' 26 Mg, 28 Si/ 29 ' 30 Si, 3 2 S/ 33 < 34 ' 26 S, 40 Ca/ 42 ' 44 > 48 Ca, 

48 Ti/ 46,47,49,50 Tii and 64 Zn/6 6,67,68,70 Zn continuously de _ 

creases toward higher metallicity. The slope changes at 
[Fe/H] ~ —2.5 and ~ —1.5 are due to the onset of 
intermediate- and low-mass AGB stars, respectively. The 
rapid change in the slope at [Fe/H] ~ — 1 is caused by SNe 
la. In Table [3] we show the mass fractions of isotopes at 
[Fe/H] = —2.6, —1.1, and 0.0, corresponding to the average 
of metal-poor SNe II, SNe 11+ AGB, and SNe II+AGB+SNe 
la, respectively, which can be compared with meteoric data 
in a future study. The results are summarized as follows. 
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Figure 17. Evolution of isoto pe ratio s again st [Fe/H] fo r the solar neighbo urhood (solid lines) with AGB yields. Observation al data 
sources include: For C, ICarretta et al.l d2000h, diamonds; ISpite et al.l (|2006h . asterisks; For Mg, lYong. Lambert, fe Ivans] d2003l), open 
and fi lled circles for 25 Mg and 26 Mg; iMelendez fe Cohenl (|2007h . open and filled squares for 25 Mg and 6 Mg; For Ti. lchavez fc Lamberj 
(2009), crosses. The solar ratios jAnders fc Grevessjl 19891 ) are shown with the solar symbols at [Fc/H] = in the upper panels. 



• Carbon — The 12 C/ 13 C ratio is ~ 4200 at [Fe/H] 
= —2.6, and decreases to 109 at [Fe/H] = —1.1 because 
of the production of 13 C from 4 - 7M© AGB stars. Then, 
because of the production of 12 C from 1 — 4M© AGB 
stars, the ratio increases until [Fe/H] ~ —0.5, and then 
decreases to be 88.6 at [Fe/H] = 0, which is consistent 
with the solar ratio of 89.9 (solar symbol, AG89) and 89.4 
jAsplund et al.ll2009l . hereafter AGSS09). The isotopic frac- 
tion of 13 C is 0.02%, 0.84%, 0.75%, and 1.03% at [Fe/H] = 
—2.6, —1.1, —0.5 and 0, respectively. The other dots (crosses 
and diamonds) show the observational data of meta l-poor 
unevolved stars (|Carretta et aUl2000l : ISpite et alj2006h . The 
low carbon isotopic values of these stars suggests that 
intermediate- mass AGB stars and/or rotating massive stars 
have contributed to galactic chemical evolution at a very 
low metallicity. Note that in the low-mass stellar models, 
non-standard extra mixing processe s are not included (e.g. 
iNollett. Busso. fc Wasserburd l2003h . These processes may 
occur during the first and asymptotic giant branches and 
result in an increase the yields of 13 C (and 14 N). While the 
effect of such processes may be minimal on the yields at 
solar metallicity, they may be substantial for Z < 0.0001, 



resulting in much smaller predicted 12 C/ 13 C ratio s from the 
AGB models than given by the yields of lKarakasI (|201Ch . 

• Nitrogen — Although 15 N destruction by AGB stars is 
neglected in our models, the 14 N/ 15 N ratio at [Fe/H] = is 
still predicted to be larger than the ratio provided by AG89 
(272). At [Fe/H] = 0, the isotopic fraction of 15 N is 0.04%, 
which is 6 — 10 times smaller than the solar ratio: 0.37% in 
AG89 and 0.23% in AGSS09. Note that the proto-solar neb- 
ula value for the 14 N/ 15 N ratio in AGSS09 is 447, a factor 
of 1.6 higher than the AG89 value. This is because AGSS09 
adopt the Jupiter nitrogen isotope value as the proto-solar 
value, noting that the AG89 ratio is the terrestrial value 
derived from air and has likely experienced isotopic frac- 
tionation and an in crease in the abundance of 15 N (see 
iMeibom etalll200l for more details). That our predicted 
nitrogen isotope value is too high is probably because the 
effect of novae, which likel y produced a substantial fra ction 
of the 15 N in the Galaxy (|Romano fc Matteucclll2003l ). are 
not included in our models. In novae, the accreted hydrogen 
is heated up to ~ 2 — 3 x 10 s K, where the CNO cycle is lim- 
ited by /3-decays rather t han the proton capture ra te of 14 N 
(the hot CNO burning. IWallace fc Wooslevl Il98ll ). There- 
fore, 13 C, 15 N, and 17 are over-produced with respect to 
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Figure 18. Evolution of isotope ratios against [Fo/H] for the solar neighbourhood with SNe II, HNe, and SNe la only (dashed lines), 
with AGB stars (solid lines), and with rotating massive stars at Z = (dotted lines). See Fig |17l for the observational data sources. 



the solar abundances. The nucleosynthesis yields depend on 
the mass of the CO and ONe white dwarfs, and the mixing 
levels between the acc reted envelope and the white dwarfs 
|jose fc Hernanj|l99ct ). The rate of novae in the Galaxy i s 
estimated to be about 30 per year (|Jose fc Hernanj|l998t ). 
but the time evolution of the rate is uncertain. 

• Oxygen — Both the 16 0/ 17 and 16 0/ ls O ratios 
rapidly decrease from ~ 190000 and ~ 14000 at [Fe/H] 
= -2.6 to 1787 and 457 at [Fe/H] = because of the metal- 
licity dependence of massive star nucleosynthesis. The iso- 
topic fractions are ( 16 0: 17 0: 18 0) = (99.99, 0.0005, 0.006), 
(99.96, 0.006, 0.037), (99.87, 0.026, 0.109), and (99.75, 0.053, 
0.194) at [Fe/H] = -2.6, -1.1, -0.5, and 0, respectively. At 
[Fe/H]=0, the 18 fraction is consistent but 17 is too large 
when compared with the solar ratios, (99.76, 0.038, 0.201) in 
AG89 and (99.76, 0.038, 0.200) in AGSS09. 17 and ls O are 
mainly produced by le O(p, 7) 17 F(/3 + ) 17 in the H-burning 
layer and 14 N(a, j) 18 F(/3 + ) 18 in the He-burning layer, re- 
spectively, and thus their abundances depend on the amount 
of the seed element. 18 O is mainly produced by He-burning 
in massive stars and is slightly destroyed in AGB stars. 17 
is over produced by AGB stars. The model that adopts only 
supernova yields (dashed line in Fig. [15)1 produces a 16 0/ 17 
ratio that is consistent with the solar ratio. Including the 
contribution from AGB stars lowers the predicted le O/ 17 



ratio to a value lower than the solar ratio. Note that for 
the majority of AGB stars in galaxies (Z 0.004), the 17 
yield is increased with the new reaction rates (§2). The oxy- 
gen isotopic ratios may be used to put constraints on the 
rates of the 17 0(p,a) 14 N and 17 0(p,7) 18 F reactions in AGB 
nucleosynthesis models. Novae also produce some 17 (and 
13 C), which would worsen the situation. 

• Neon — The 20 Ne/ 21 Ne and 20 Ne/ 22 Ne ratios show a 
similar decrease as oxygen. With AGB yields, the isotopic 
fractions at [Fe/H]= are consistent with the solar ratio: 
( 20 Ne: 21 Ne: 22 Ne) = (96.16, 0.23, 3.61) in the model, (92.99, 
0.23, 6.78) in AG89, and (92.94, 0.22, 6.83) in AGSS09. The 
only data for comparison is meteoritic. 

• Magnesium — The 24 Mg/ 25 Mg and 24 Mg/ 26 Mg ratios 
also show a rapid decrease from 247 and 271 at [Fe/H] 
= -2.6 to 8.54 and 9.22 at [Fe/H] = because the produc- 
tion of the minor isotopes increases in metal-rich supernovae. 
The 24 Mg/ 25 Mg and 24 Mg/ 26 Mg r atios are a bit larger than 
the observations of stars (dot s, lYong. Lambert, fc Ivansl 
120031 ; iMelendez fc Cohen! 120071 ) and the solar ratios (7.92 
and 7.19 in AG89), which suggests that AGB stars (or 
Wolf-Rayet stars) need to contribute more at all metal- 
licities. The isotopic fractions are ( 24 Mg: 25 Mg: 26 Mg) = 
(99.28, 0.39, 0.34), (95.60, 2.30, 2.11), (88.09, 5.77, 5.13), 
and (82.46, 9.28, 8.26) at [Fe/H] = -2.6,-1.1,-0.5, and 
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Figure 19. Evolution of isotope ratios against [Fe/H] for the solar neighbourhood (solid lines), halo (short-dashed lines), bulge (long- 
dashed lines), and thick disk (dot-dashed lines) with AGB yields. See Fig |17l for the observational data sources. 



0, respectively. The isotopic fractions of solar ratios are 
(79.03, 9.97, 10.99) in AG89 and (78.99, 10.00, 11.01) in 
AGSS09. In both observations, 26 Mg is more abundant 
than 25 Mg, which is the opposite both for the AGB stars 
and supernovae yields. In supernovae, there is a metallic- 
ity dependence on the predicted Mg isotopic ratios. 24 Mg 
is mainly produced by the reaction 12 C+ 12 C and is a pri- 
mary isotope. On the other hand, 25 Mg and 26 Mg are sec- 
ondary isotopes and are significantly produced by higher 
metallicity stars. They are mainly produced by the follow- 
ing reactions: 24 Mg(p, 7 ) 25 Al(^+) 25 Mg(p, 7) 26 Al(/3+) 26 Mg 
and thus their abundances depend on the amount of 
the seed element, 24 Mg. In AGB stars, there is no such 
metallicity dependence. The neutron-rich Mg isotopes are 
produced by cv-cap tures onto 22 Ne via 22 N e(a, w) 25 Mg 
and 22 Ne(«,7) 26 Mg (|Karakas fc Lattanziolliool ). The AGB 
yields of the neutron-rich Mg isotopes depends on the 
amount of 22 Ne in the He-burning region, and this can have a 
primary component in intermediate-mass AGB stars, where 
the 22 Ne is produced from primary nitrogen made by HBB. 

• Si, S, Ca, Ti, and Zn — For heavier elements most of the 
ratios at [Fe/H] = are roughly consistent with the solar ra- 

29 48 47 

tios, but the predictions for Si, Ca, Ti are smaller than 
the solar ratios. These may require the update of the mixing 
treatment and reaction rates in the supernovae calculations. 



For Zn, the offsets from the solar ratios are possibly caused 
by the under-production of 64 Zn. Although neutron-rich iso- 
topes of Zn could be produced by neutron-capture processes, 
64 Zn is mostly produced by the higher energy and entropy 
experienced during supernova explosions. In fact, if we were 
to set a higher fraction of hypernovae, the isotopic ratios 
become closer to the solar ratios. 



The model dependence on the evolution of the isotope 
ratios is shown in Figure 1181 With the contribution from 
AGB stars (solid lines), the 12 C/ 13 C and 20 Ne/ 22 Ne ratios 
become consistent with the solar ratio, and 24 Mg/ 25 ' 26 Mg 
becomes slightly closer to the observations. However, as 
mentioned above, the le O/ lr O ratio is consistent with the 
model withou t an AGB contribution (dashed line s). This is 
different from lTimmes. Wooslev. fc Weaverl (|l995l )'s results, 
where the 16 0/ lr O ratio was consistent with the solar ra- 
tio without AGB yields and 18 O is over-produced by super- 
novae. This is due to our use of updated supernova yields. 
For heavier elements, there is only a ~ 0.01 dex difference 
around [Fe/H] ~ for the models with and without AGB 
yields. F or 24 Mg/ 25 ' 26 Mg, th e contribution of AGB stars are 
larger in iFenner et ail (2003) than in our models, which may 
be due to differences in the adopted IMF (§2.3). 
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If we include rotating massive stars (dotted lines), the 
12 C/ 13 C ratio becomes much closer to the observatio ns of 
metal-poor stars as shown by IChiappini et al . (200|), but 
16 0/ 17 cannot be improved when comparing with the so- 
lar ratio. Note that the N and Mg yields of rotating massive 
stars are not included. The small difference among the mod- 
els with and without rotating massive stars is caused by the 
difference in the IMF (M Ul 2 = 120Mq is adopted in the case 
with rotating massive stars instead of Af Ui 2 = 50Mq). 

The dependence on the star formation history is shown 
in Figure Ql] For 12 C/ 13 C, the first decrease is due to 13 C 
production from intermediate-mass AGB stars, the next in- 
crease is due to the 12 C production from low- mass AGB 
stars, and the following second decrease is due to SNe la. 
These modulations appear at a higher metallicity in the 
bulge (long-dashed lines) and thick disk (dot-dashed lines) 
than in the solar neighbourhood (solid lines). On the other 
hand, these evolutionary changes appear at lower metallicity 
in the halo (short-dashed lines). In general, the ratios be- 
tween the major and minor isotopes such as 24 Mg/ 25,26 Mg 
are smaller in the bulge and thick disk, and are larger in the 
halo because of the metallicity effect of supernovae. How- 
ever, the 16 0/ 17 ratio in the halo is low due to the pro- 
duction of 17 O from low-mass AGB stars, as also seen in the 
high [(C, F)/Fe] abundances. Therefore, the isotopic ratios 
can be used as a tool to pick out the stars that form in a 
system with a low chemical enrichment efficiency. This may 
be possible in our halo, but more likely in small satellite 
galaxies that were accreted onto our Milky Way Galaxy. In 
Table|31 the mass fractions of isotopes at [Fe/H] = —0.5 are 
provided for the solar neighbourhood, halo, bulge, and thick 
disk models. 



4 CONCLUSIONS 

We have presented evolution of isotope ratios of elemen- 
tal abundances (from C to Zn) in the solar neighbourhood, 
bulge, halo, and thick disk, using chemical evolution mod- 
els with the updated yields of AGB stars and core-collapse 
supernovae. Although the ejected mass of heavy elements is 
much smaller for AGB stars than for supernovae, the AGB 
yields of C, N, and F are comparable to those of super- 
novae. When examining the AGB yields, the isotopes of 12 C 
and 19 F are mainly produced in low-mass stars (1 — 4Mq), 
whereas 13 C, 14 N, 22 Ne, 25 Mg, and 26 Mg are mainly pro- 
duced in intermediate- mass (4 — 7Mq) stars. By including 
AGB yields, the predicted ratios of [C/Fe], [N/Fe], 12 C/ 13 C, 
and 20 Ne/ 22 Ne at [Fe/H] = are much improved compared 
with the observed ratios in the solar neighbourhood. For su- 
pernovae, minor isotope production increases with increas- 
ing metallicity. Therefore, the [(C, F)/Fe] and isotope ratios 
can be used as a cosmic clock along with the [a/Fe] ratio. 

Because of the effect of the progenitor mass and metal- 
licity, the evolutionary history of elements varies for environ- 
ments with different star formation histories. In the bulge 
and thick disk, the star formation timescale is shorter than in 
the solar neighbourhood, and [a/Fe] is higher and [Mn/Fe] 
is lower because of a lower contribution from SNe la. In 
contrast, the [(Na, Al, P, CI, K, Sc, Cu, Zn)/Fe] ratios are 
higher because of the effect of metallicity. In other words, 
from these elemental abundance ratios, it is possible to se- 



Table 2. The nucleosynthesis yields of core-collapse supernovae 
updated from Kobayashi ct al. (2006) in the ejecta in Mq. 
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4.13E-06 
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7.68E-06 


48 Ti 


2.49E-04 
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1.01E-05 
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1.04E-06 


1.51E-05 


1.32E-05 


50 V 


1.86E-08 


9.94E-08 


2.07E-07 


Sly 


1.57E-05 


1.20E-05 


1.80E-05 


50 Cr 


5.15E-05 


5.04E-05 


8.15E-05 


52 Cr 


3.89E-03 


3.50E-03 


1.26E-03 


53 Cr 


2.40E-04 


9.95E-05 


1.06E-04 


54 Cr 


2.75E-06 


3.00E-05 


2.70E-05 
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0.02 
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M 


18 


25 


25 


E 


1 


1 


10 
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54 Fe 
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57 Fe 
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S 41 F 04 

O .t: J. Ill Urfc 


59 Co 


3 58F-05 


4.01E-04 
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58 Ni 
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63 Cu 
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9 90F 04 
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7.19E-06 


5.65E-05 


2.88E-04 


66 Zn 


2.35E-05 


3.34E-04 


3.47E-04 


67 Zn 


3.16E-06 


8.10E-05 


7.26E-05 


68 Zn 


2.37E-05 


4.63E-04 


4.24E-04 


™Zn 


4.95E-07 


1.14E-05 


1.94E-05 


69 Ga 


2.30E-06 


5.77E-05 


5.19E-05 


71 Ga 


2.53E-06 


3.36E-05 


3.61E-05 


70 Gc 


3.08E-06 


7.28E-05 


6.25E-05 


72 Ge 


5.07E-06 


1.19E-04 


1.07E-04 


73 Ge 


6.71E-07 


2.73E-05 


2.32E-05 


74 Ge 


9.69E-06 


3.34E-04 


2.87E-04 



lect metal-rich stars ([Fe/H] > — 1) that formed in the early 
Universe as a result of rapid star formation. 

In the halo, the chemical enrichment timescale is longer 
than in the solar neighbourhood, and [(C, F)/Fe] and 
12 C/ 13 C are higher because of a stronger contribution from 
low- mass AGB stars. The [a/Fe] and [Mn/Fe] relations are 
the same as in the solar neighbourhood, but [(Na, Al, P, CI, 
K, Sc, Cu, Zn)/Fe] are lower and other isotope ratios such 
as 16 0/ 18 and 24 Mg/ 25 ' 26 Mg are larger because of the low 
metallicity. Not only from elemental abundance ratios, but 
also from isotopic ratios, it is possible to select the stars 
that formed in a system with a low chemical enrichment 
efficiency. 

Isotopic ratios provide useful information to improve 
theoretical predictions in terms of reaction rates, nucleosyn- 
thesis, and the modelling of convective mixing. The under- 
production of N suggests the contribution from novae, the 
over-production of 17 requires the updating of reaction 
rates that feed into AGB yields, and the under-production 
of 29 Si, 48 Ca, 47 Ti and possibly 64 Zn may require the up- 
dating of supernova models. The observational constraints 
at low-metallicity are particularly important. However, it is 
not easy to estimate isotopic ratios from observations of stel- 
lar spectra. To detect the small line shifts between isotopes, 
very high quality data is required. The next generation of 
large telescopes (e.g., the Giant Magellan Telescope and the 
European Extremely Large Telescope) will be constrained 
by the same problems but will be able to measure isotopic 
ratios out to larger distances. It will then be possible to 
study the evolution of isotopic ratios in the local neighbour- 
hood of our Milky Way Galaxy including in dwarf spheroidal 
galaxies. These data will help to put important constraints 
on the star formation and chemical enrichment histories of 
such systems, and which may lead to the verification of the 
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hierarchical clustering in CDM cosmology using chemody- 
namical simulations. 
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Table 3. The mass fractions of isotopes for the models of the solar neighbhorhood, halo, bulge, and thick disk. 
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[Fe/H] -2.6 -1.1 -0.5 -0.5 -0.5 -0.5 
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6 Li 
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7 Li 
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7.82E- 


■03 


6.09E-03 


17 o 


4.23E- 


■10 


1.41E-07 
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09 
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40 K 
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10 
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■06 
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■05 
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